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Abstract 

We derive the equation of state of symmetric nuclear matter in a relativistic 
■ theory with a and lo exchange. We take a chiral version of this model which insures 

all the chiral constraints. Going beyond the mean field approach we introduce the 
^Sj ' effects of the pion loops. For the parameters of the model, in order to fix those linked 

^ ■ to pion exchange, we exploit the most recent information on the short-range part of 

O . the spin-isospin interaction. For those linked to the scalar meson exchange we make 

I use of an analysis of lattice results on the nucleon mass evolution with the quark 

^T) ' mass. With these inputs we are able reach a correct description of the saturation 

properties. From the equation of state of symmetric nuclear matter we alsoderive 
the density dependence of the quark condensate and of the QCD susceptibilities. 
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1 Introduction 

Relativistic theories of nuclear matter are needed for the prediction of the equation of 
state of nuclear matter at large densities, or of other quantities such as the QCD quark 
^ , condensate or the QCD susceptibilities. They find a natural framework in the relativistic 
^ I a and u exchange model [I]. One important aspect is its chiral nature, especially if it 
^ ■ has to be used also for the investigation of QCD related quantities such as the quark 
condensate. The approach that we follow in this work for symmetric nuclear matter is 
^ . based on the chiral version of the a and u model that we have proposed in ref. [2]. 

Here the sigma model is formulated in a non-linear version but with the presence of a 
chiral singlet scalar field s. The nuclear binding is provided by this chiral invariant field 
and in this way all chiral constraints are automatically satisfied. However it is known 
that in such chiral theories attractive tadpole diagrams destroy saturation [31 H]. Hence 
in order to achieve saturation we complement our model, as in our previous work [5], 
by the introduction of the scalar nucleonic response first introduced in the quark-meson 
coupling model, QMC [6j. In our previous investigation [5], we have applied this model 
in the Hartree approximation to derive the equation of state. From this we have deduced 
the density evolution of the quark condensate and of the QCD susceptibilities. This 
procedure insures the compatibility with the saturation properties of nuclear matter. Our 
conclusions were that the deviations of the condensate evolution with density from the 
linear behavior are mild and that there exists a convergence effect between the scalar and 
pseudoscalar susceptibilities, a conclusion already reached in ref. [TJ [8] . 

The pion which does not contribute at the mean field level was ignored in this work. 
However the exchange of the scalar s field between nucleons is not the only source of 
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attraction. Beside the exchange of a scalar meson, the Van Der Waals type of forces, i.e., 
two-pion exchange with one or two intermediate Delta excitations, also contribute to the 
middle range attraction of the NN potential and to the nuclear binding. Iterated pion 
exchange is also a source of attraction. It is even envisaged in a chiral theory that it could 
describe the bulk of the nuclear binding [9J • It is therefore natural to extend our previous 
descriptions restricted to the Hartree scheme so as to incorporate pion loops, which is the 
aim of the present article. When pion exchange forces are discussed one has also to intro- 
duce the other components of the spin-isospin force, namely the short-range contribution 
embodied in the Landau Migdal parameters. We will also introduce rho exchange, which 
acts in the spin transverse channel. We will show that the role of both the short-range 
terms and the rho in the correlation energy is crucial. For what concerns the values of 
the Landau-Migdal parameters, years of experimental and theoretical investigations have 
sharpened our knowledge of these quantities. We will exploit here the most recent results, 
for which a survey is given in a recent article by Ichimura et al. [10] . We feel timely to 
revisit certain problems of nuclear physics and combine in these studies the experience ac- 
cumulated in the pionic type of physics including chiral aspects and some results of QCD 
to have simultaneously a description of the nuclear binding and the evolution of chiral 
symmetry restoration, i.e., of the quark scalar density, and of the QCD susceptibilities. 

In order to achieve this we will incorporate in a consistent way, on top of this mean 
field, the pion loop contribution to the energy, namely the pion Fock term and the pionic 
correlation energy, taking into account short-range correlation effects, within a RPA ring 
approximation [11] . In this way, by taking the successive derivatives of the grand potential 
with respect to the pion mass, we will obtain the full contribution to the evolution of 
the chiral condensate and the scalar susceptibility. In a recent work [12] we obtained 
the pionic piece of the scalar susceptibility by taking the derivative with respect to the 
quark mass (in practice to the pion mass) of the pion scalar density ($^), which is a 
component of the condensate of the medium. We showed that the medium effects affecting 
the scalar susceptibility are precisely the same that the ones governing the reshaping 
of the strength in the low invariant mass region in two-pion production experiments 
[131 HH [12]. One purpose of the present paper is to provide a unified description of 
two quantities, i.e. the effect of the nuclear scalar field and of the pion loops in a way 
compatible with saturation properties of nuclear matter. The pionic contributions depend 
on the pionic polarization propagator i.e. the p.h. (particle-hole) and A.h. (Delta-hole) 
Lindhardt functions. Although we start from a relativistic mean field theory we will 
evaluate these Lindhardt functions in a non-relativistic framework for convenience since 
it is based on a well established phenomenology. It is interesting to point out that this 
approach where pion loops are incorporated on top of background scalar and vector mean- 
fields has some similarities with the approach of reference |16j based on a relativistic 
density functional. There are however important differences that will be discussed in 
section 4. 

Our article is organized as follows. Section 2 is devoted to the loop energy. After a 
brief summary of the results obtained in the Hartree scheme pion loops are introduced 
together with the other short-range components of the spin-isospin interaction. Section 3 
discusses the QCD quantities, the quark condensate and the scalar susceptibility. Finally 
in section 4 numerical results are given and the results discussed. 
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2 Pion loop energy 



2.1 Summary of the mean field results without pion loop 

We start with the mean field Hamiltonian used in ref. ||5j supplemented by the free pion 
Hamiltonian : 

where s is the chiral invariant field associated with the radius 5 = /tt + s of the chiral 
circle. V{s) = Vo(s) — cS is the vacuum potential which can be split into Vo(s) = 
(A/4)((/^ + s)^ — v^) , responsible for spontaneous chiral symmetry breaking, and the 
explicit symmetry breaking piece, ~cS, where c = f-^fn^- As usual in relativistic mean 
field theories we have added a coupling to an omega field ujq . 

At the mean field level the energy density for symmetric nuclear matteris given by : 



and = ^Jp^ + M*^{s) is the energy of an effective nucleon with the effective mass 

M]^(s). The effective mass is sensitive to the effect of the nucleon susceptibility, i^ns^ 
which embeds the influence of the internal nucleon structure : 

M^(S) =Mn {l + J^ + \ KNS s\ (3) 

The expectation value, = /^r + s, of the S field plays the role of a chiral order parameter. 
It is obtained by minimizing the energy density : de/ds = g*s Ps + ^'(^) = with the 
following expressions for the scalar density, ps-, and the scalar coupling constant g*^ : 



The quantity gs = Mj^/f^^ is the scalar coupling constant of the model in the vacuum. 
Notice that the density dependence of g^ arises entirely from the susceptibility term. 
Since the mean scalar field is negative and the sign of k^s is positive, g^ is a decreasing 
function of the density. The in-medium sigma mass is obtained as the second derivative 
of the energy density with respect to the order parameter : m*^ = d'^e/ds^ = V"{s) + 
d{9*sPs)/ds. 

In the free pion part of the Hamiltonian of eq. [H the operator b]^ creates a free pion 
state with momentum q and isospin index j and we we have Uq = ^/q"^ + where m^r 
is the pion mass. In fact in a more complete treatment including s-wave pion-nucleon 
coupling, according to ref. [2], the pion mass should be replaced by an effective mass 
m*^ = + s/fir). As in ref. [T7j we do not consider this effect since it is almost 

completely compensated by other contributions [18] and the pion mass remains stable in 
the nuclear medium [19]. 
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2.2 Inclusion of the pion-nucleon coupling 



We now include the usual derivative pion-nucleon interaction Hamiltonian : 

H^NN = - j dr^NYl^^A^N ^- J dr^N^(T.V$.rN (4) 

where the second form corresponds to the non relativistic limit of pure p-wave nature. It 
can be rewritten in a standard second-quantized form as : 

1/2 



(1 \ 



where we have included phenomenologically a dipole nNN form factor f (q) and V is 
the volume of the normalization box. Using Green's function techniques, the expectation 
value on the nuclear ground state |0) of this interacting Hamiltonian can be obtained 
according to [11] : 
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{H^^^) = 3V — ^ UUuJ, q). (6) 

Here H2,(ci;,q) is the full spin-isospin polarization propagator defined by : 

H^(^,q) = ^^°°d(t-t')eM*-*')(z!)(o|T(L](q,t),L,(q,t'))|0) 

1 y Kn|L,(q)|0)p |(n|Lt(q)|0)p ^ 
V ^ uj - En + ir] u + En - ir] 

n 

We also include the Delta-hole states. In this case the operator cr.qr, is replaced by 
{.Q-kN/^I 9-kNn) S.q Tj where the Sj (Tj) are the standard spin (isospin) transition operators 
between spin (isospin) 1/2 and 3/2 states. For the ratio of coupling constants, we take 
the quark model value : -Rata = {q-kNh/ Q-kNn) = \/72/25. 



2.3 Inclusion of short-range correlation and rho meson exchange 

The expression ([6]) written above shows that the interaction energy of the pion cloud 
with matter is obtained by integrating the product of the full longitudinal spin-isospin 
polarization propagator {i.e. the nuclear response to a pion-like excitation) with the pion 
exchange potential. However, due to the correlated nature of the medium, the contact 
piece of the pion exchange potential should be treated with caution. In the part which 
concerns medium effects, this piece should be removed. This is achieved as usual by 
introducing in the spin-isospin channel a short-range two-body potential governed by 
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Landau-Migdal g' parameters. The Nh — Nh piece of the corresponding Hamiltonian 
writes : 



Hsr = ^ J]ci4c5C^(a/3|K,(l,2)|75) 

K.(l,2) = 1^ (||yt;2(q)^U'T(l).^(2)ni).n2)e-''>-(^-^^^ (8) 

It is straightforward to extend this expression in order to include the Nh — Ah and the 
Ah — Ah pieces. For this we replace g'j^^ by RnaQna -^a^aA'aa where -Rata defined 
previously is the coupling constant ratio. In the numerical calculation we will use the 
following set of values suggested in ref. [10], which takes into account the most recent 
data analysis : g'^j^ = 0.7, (7^^ = 0.3, g'/^^ = 0.5, with clear deviations from universality 
as advocated in ref. [10]. The spin-spin operator can be decomposed into a longitudinal 
piece and a transverse one. For this reason we also introduce the full transverse spin- 
isospin polarization propagator : 

nT(a;,q) = ^|_^°°rf(t-Oe^^(*-*')^(0|T(Tt(q,^^ 

1 |(n|T,(q)|0)p |(r^|Tt(q)|0)p 

2V ^ uj - En + W uj + En - iri ^ ' 

n 

where Tj is obtained from Lj by replacing cr.q by cr x q. The expectation value of the 
short-range potential is : 

{Hsr) = \v J^^'-^J-^^g'[{U,{uj,ci) + 2UT{uj,ci))] - ^ (a/.|K.(l, 2)|M 

(10) 

where the second term is introduced in order to remove the contribution for free nucleons. 
In the following, when this second term term is not explicitely written, this substraction 
is done implicitely. The first term g' {IlL{uJ,q) + 2nT(ci;,q)) is written in a schematic 
form. The precise form involving the three Landau-Migdal parameters, g'NNy 9NAy dAA^ 
will be given below in subsection 2.4. 

Finally, so as to obtain a more realistic description of the transverse spin-isospin 
channel we also introduce, beside the short-range piece, the p meson through a pNN 
Hamiltonian. In the non relativistic limit it has a structure similar to the irNN one, once 
the cr.q coupling is replaced by a cr x q coupling. Its expectation value is : 

^•^idu f dq 2Cp 



{ff,™) = 3V' / _ p^;3^n.(..q), (11) 



with fig = y/g2 -|- m'j^. We take Cp = 2 which corresponds to the strong rho coupling [20j. 
2.4 Calculation of the pion loop energy 

In order to calculate the pion loop energy we use the well-known charging formula method. 
For this purpose we introduce an auxiliary Hamiltonian which depends on one strength 
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parameter A : 

H{X) = Hq + X {Ht,nn + HpNN) + A ■ (12) 

H{X) coincides with the original Hamiltonian for A = 1. The charging formula allows the 
evaluation of the modification of the ground state energy with respect to the ground state 
energy Eq of the Hq Hamiltonian : 

E,,,,-E-Eo = Ve,..,--V ^ j (^ X T 

([rz.(cc',q)ni(cu,q;A)] + 2 [^^(cu, q) 0^(0;, q; A)] ) . (13) 

The dependence on A of the full polarization propagators can be obtained by systemati- 
cally multiplying the coupling qa/ fn by A. For instance there appears a net A^ factor in 
front of each first order Lindhardt function appearing in the RPA expression of the full 
polarization propagators. The quantities Vl,t(i^, q) are the full energy dependent residual 
interactions in the longitudinal and transverse spin-isospin channels. In the Nh — Nh 
sector, they write : 

Vljv7v(^, q) = q'nn + —2 2^ Vtnn{^, q) = q'nn + Cp — (14) 

with identical expressions in the Nh — Ah and the Ah — Ah channels once the relevant g' 
parameters have been incorporated. Again the expression (fT3|l written above is schematic 
and the full formula calculated in the RPA ring approximation is given just below in eqs. 
[WT71 The A integration can be done analytically to obtain the loop energy density eioop '■ 

Sloop = Sl + ex with (15) 
^L = lj [-ln{l- VlnnII% - VlnaH'a - (V^jva - Vlnn Vlaa)K^V)] (16) 

^T = sJ '-^^ [-ln{l- VtnnH'n - VtnaHI - (VSna - Vtnn Vtaa)II%III)] (17) 

is the N.h. Lindhardt function which is calculated non relativistically. With the 
single particle energy, Sp = /2M* where M* is the s dependent nucleon effective mass 
calculated at the mean-field level, 11^ has the familiar form : 

n(o.,q)=4f^V.^(q) ^ fQ{pF-p)Q{\v + ^\-PF) 



2/^7 J (27r)3 V UJ - eivpq + irj 

Q(Pf -p)B(|p + q| -Pf) 
+ e^Ypq - iv 

eATpq = ep+q — ep is the energy of the p.h. excitation. The corresponding A.h. Lindhardt 
function, 11^ is also calculated non relativistically. Introducing also an effective A mass, 
we make the natural assumption that the in-medium shift of the nucleon and delta masses 
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in presence of the nuclear scalar field are identical. The A.h. energy is thus : CApq = 
eA,p+q — Cp with eA,p = Ma — +p'^ /2M^ and 11^ writes : 

n^(-,q) = ^ fif ^-^V^'^k) / ^e(p.-p)f . 

9 \2f^ J J {2TTy \u - CApq + tr] u + CApq/ 

(19) 

This loop energy can be split in a piece of mean-field nature (containing the Fock term) 
and a genuine correlation energy. For instance we have for the longitudinal piece : 

^MF I ^Corr 

2 / (27r)^ [VLATArll^ + VlataHa] 

I j [-H^- VlnnK - VlnahI - (\/^V - Vlnn VLAA)n^n^) 

-VlnnII% - VlnaT^I] ■ (20) 

By performing the energy integration with a Wick rotation, e^^ can be expressed as 
follows : 

^r = /(^0(Pi^-p)s.(p) + ef-^ (21) 

S^(p) is the pionic contribution to the nucleon self-energy : 

S-(P) = -o hrr / TTT^q ^ 7 + — 7 • (22) 

2 V2/vr/ i (27r)'^ V^g^^g + eAfpq 9 ujqUJg + eApqJ 

In principle medium effects are present in this quantity since the N.h. and A.h. energies 
depend on the effective nucleon and delta masses but in practice we ignore these effects. 
The momentum dependence of the nucleon self-energy is linked to the treatment of the 
vacuum nucleon which is not the purpose of this paper. Hence we ignore it. At zero 
momentum the self-energy = S^(p = 0) is similar to the model calculation one used 
to fit lattice data by Thomas et al [21], although we do not make a static approximation. 
In our calculation we will choose a dipole form factor, with a cutoff A = 0.98 GeV, such 

(it) 

that the resulting contribution to the free nucleon sigma term, cr]^ , the expression of 
which being given in subsection 3.1, is cx^'' = 21.5 Me V^, in agreement with previous 
determinations [221 [23] . The corresponding pion cloud self-energy is T,.^ = —420 MeV. 
The Fock term is : 

{23) 

Notice that this Fock term includes a retardation effect through the presence of the 
particle- hole energy CATpq. In practice the effect is numerically small and we take it as its 
average Fermi sea value : eA^pq = e^g = q'^/2M]y, ignoring the in-medium modification 
of the nucleon mass. In that way the Fock energy does not depend on s, a feature which 
simplifies the calculation of the scalar susceptibility. Also notice that, at variance with 
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other works, see e.g. [23], we keep the form factor in the Fock term. The Fock term 
relative to the transverse channel is obtained in a similar way : 

(24) 

Finally the correlation energy is also calculated using a Wick rotation for the u energy 
integration. The numerical results will be given in section 4. 



3 The chiral condensate and the scalar susceptibihty 
3.1 Direct calculation 

In this first approach the in-medium quark condensate is evaluated directly from the 
expectation value of chiral symmetry breaking piece of the Hamiltonian : 

(H^sb) = V2m{qq) ^ V (^^ml^^ - cS^ ^ {qq) ^ {qq)vac - Jjl + f) 

'(25) 

where the GOR relation has been used to obtain the second equation. The pion scalar 
density can be directly calculated from the in-medium pion propagator [25l [12] : 



{^') = 3 1 ^ [D^ - Do^] {uj, q) = 3 1 ^ [dm-^, q) n°(^, q)] • (26) 

Dj^ = (cij^— tu^— q^n°)^^ is the full pion propagator, q^n° is the irreducible pion self-energy 
with : 

j^O ^ + + i'^dNA ~ 9'nN ~ 9'aa)^N^A ^27) 

1 - + 9'aaK + «A - g'NNd'AAm^i) ■ 

It is possible to show that the formal expression of the pion scalar density can be obtained 
as the derivative with respect to of the pionic energy density El given by eq. 



-^ = d^^ or equivalents -{mU.jjj = -g^ (28) 

which is nothing but the Feynman-Hellman theorem applied to the pion loop contribution. 
As in the previous section we decompose into a free nucleon contribution, a Fock 
contribution and a correlation piece, with a similar decomposition for ($^) : 

= p ^ + M!!! + M!I ^ p ^ + + (29) 

2 (9m^ dm'l dm'l 2 2 

The pionic contribution to the pion-nucleon sigma term is : 



W _ 3 / V 2 f q^v^{q) 



- 2 [2fJ J (2vr)3 2^,2 



11 1 

^_UJq UJq + eg {Ujq + e,)^ 



9 \ujq UJq + eAq (cUq + eAqf 
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(30) 



The Fock term contribution has the following explicit form : 
2 ~~2V2aJ ^ Jo (27r)3 2cu2(u;, + e,) 



Finally the correlation piece is obtained by removing the leading order piece in eq. [26] : 

~2~ =2 J 72^ " " ^1 "l^- ^^^^ 

In the next subsection we will evaluate the condensate from the equation of state (grand 
potential) and the Feynman-Hellman theorem. It leads essentially to the same result with 
a modification coming from the fact that we will have to face the problem of the Lorentz 
nature of the in-medium effects of the pion loop energy. 



3q_ 

ApF 



+ 



1 

— + 



(31) 



3.2 The quark condensate from the equation of state 

The pion cloud contribution has been calculated within the standard non relativistic 
framework. One reason is the basis of the pion-nucleus phenomenology. However some 
care should be taken when this pion cloud is included on top of the relativistic mean 
field used to describe the coupling of the scalar field to nucleons. For free nucleons the 
pionic self-energy has to be of scalar nature and thus it adds to the bare nucleon mass: 
Mtv = Mo + Stt. For what concerns the in- medium contribution to the energy we do not 
know if it is of scalar or vector nature or a mixture of both. Here we assume its vector 
nature and add it to the energy density. This assumption affects contributions arising 
from the pionic interaction between nucleons where it amounts to distinguishing between 
scalar and vector densities. In practice the effect of this difference is numerically quite 
small. For the free nucleon part instead this distinction has a non negligible effect. The 
energy density thus writes : 



where e^^ = e^°'^^ + e^°'^^ + e^""^^ + Ej-""^^ summarizes the loop energy density from the 
Fock term and correlations (pion + rho + short-range) which can be decomposed into a 
longitudinal piece (pion -|- short-range) and a transverse piece (rho + short-range). The 
mean field s is obtained by minimizing the energy density before the incorporation of 
the loop energy or, said differently, ignoring the s dependance of e^^ . This is actually 
true for the Fock term and this approximation has a small numerical incidence on the 
evolution of the condensate with inclusion of the correlation energy. It also follows from 
the independence on s of e^*" that the in-medium sigma mass remains the same as in the 
mean-field case : 



< = Q-2= ^"^) + Ps + 9*s-^- (34) 
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As discussed in ref. [5] the scalar susceptibility term (second term in n^s on the r.h.s of 
eq. \3^ counterbalances the effect of the first term, the in-medium chiral dropping of the 
sigma mass. The last term actually corresponds to the nuclear response associated with 
NN excitation. In practice it is small and it can be omitted [8]. 

We will derive the in-medium chiral condensate and the QCD scalar susceptibility 
from the equation of state since they are related to the first and second derivatives of the 
grand potential with respect to the quark mass m at constant chemical potential /i. The 
baryonic chemical potential is obtained as : 

2 f)^FC 



/. = | = i?; + ^P+^ with E*,= ^pl + M*^rs)- (35) 

In order to obtain this result we have taken into account only the explicit density depen- 
dence on the density p, i.e., keeping s constant since de/ds = 0. One deduces that the 
baryonic density is controlled by the chemical potential according to : 

p = j Wf^K!'-''^-^/-^) ^''^ 

while the scalar density writes : 

The grand potential, which is obtained through a Legendre transform, can be written in 
the following form : 

uj{fi) = e - fip 



2 f) FC 

+ Vis)-^p^ + e^^ -p^. (38) 
2 op 

For the derivation of the condensate and of the susceptibility we have to specify the sym- 
metry breaking parameter, which in QCD is the quark mass. In the context of this model 
it is the quantity c = /,r^^ which enters the symmetry breaking piece of the potential. In 
the application of the Feynman-Hellman theorem we use the explicit expression of dc/ dm 
given by the model to leading order in the quark mass m, i.e., dc/dm = —2{qq)^ac/fn- 
As previously, in the calculation of the derivative we only keep the explicit dependence 
on c : 

2\dm)^ 2dm\dc)^'^ U \ ^ \ dc ) g dc) 

S ml ml del^ 



fj, c dml c dml 

imUc ( 1 + ;^ - ^) • (39) 
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In the last line we have grouped the pionic contribution in a pion scalar density given by : 



ps 



a 



N 



+ 



2\Fock 



+ 



2\Corr 



(40) 



Notice that the pion scalar density deviates from the fully non relativistic one of eq. [26] 
since the leading order term behaves like ps in place of p. To leading order in density we 
recover, as expected, the well-known result : 



(gg) = {qq)^ 



1 



ctn p 
mi 



m 



Ugs^. (41) 



3.3 The scalar susceptibility 

The in-medium scalar susceptibility is obtained as a derivative at fixed p of the chiral 
condensate given in eq. [39l The dependence of the pion scalar density on the symmetry 
breaking parameter (c or m^) receives two contributions : one comes from the explicit 
dependence on c and the other from the implicit dependence since, at fixed p, the density 
p depends on c according to eq. [361 Hence the scalar susceptibility takes the form : 



Xs 



( 9{qq) \ 1 / dc 

\ dm J 2 \ dm 



m 



-2 



2 

vac 



J IT 



fdS_ 



1 



1 



dc 



dp 



dp 
dc 



m. 



d{^^) 
2/^c V dml 



nuclear , pionloop 

Xs Xs 



(42) 



We have divided xs in two components. The first one which is the sum of the first two 
terms contains the coupling of the scalar quark density fluctuations to the nuclear p.h. 
excitations. For this reason we denote it ■y/^^'^^^'^'^ ; 



nuclear 

Xs 



-2 



2 

vac 



J TT 



dS\ 
~d~c ) 



dp 
dc 



(43) 



The second one, ■y/P^°^^°°P ^ of purely pionic nature, comes from the explicit dependence 
on of the pion scalar density ($^). The latter can be written as : 



pionloop 

Xs 



'vac 



1 9($ ) nion—Fock . pion—Corr 

- Ps Xns + Xs + Xs 



n 2/2 dml 



(44) 



The pionic contribution to the nucleon scalar susceptibility Xns obtained from the 
derivitative of the pion cloud contribution to the sigma commutator (eq. [30l) : 



(tt) 

Xns 



2 

vac 



J w 



3 
2 



9a 
2U 



4 -Rata 



dq q f (q 

(2^ 
3 



{Ujq + tq) UJq {ujq + tg^ {ujq + e,)^ 



9 \UJ^ {Ujg + CAg) 



+ 



UJg {uJg + e^gY 



{Ug + e^gf 



(45) 
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The Fock term contribution is obtained from eq. 



pion—Fock n \ii/vac - i \ „ ; — i t. - i i ""^ i 



' ' + US) 



{Ujq + e,) UJq {Ujq + tqY {ujq + tq)^ 

Finally the correlation piece is obtained from eq. 



xr-"'- = 2 ^ ^ / 1^ [( ^'^0. + D^Dl ) n° - 2 n°] {u, q). (47) 



In the evaluation of the nuclear contribution, -i^g^^ear^ derivative {dS /dc)^ is obtained 
by taking the derivative of the minimization equation with respect to the parameter c : 



m 



*2 



9s 



M* 

l^no(o) 



(4J 



with the auxiliary quantities : 
2/. 9p 



+ 



M 



N 



+ 



9^5 



(49) 



and with no(0) = —2M'^'ppj'n^^ which is the non-relativistic free Fermi gas particle-hole 
polarization propagator in the Hartree scheme, at zero energy in the limit of vanishing 
momentum. The derivative of the baryonic density is obtained by taking the derivative 
with respect to c of eq. [36], with the result : 



dp 
dc 



9s 



/as 



E 



(-) +^^^'^Mno(o)(i- wno(o)) 



-1 



(50) 



N 



It follows that T^g^cZear written in the following form which displays the propagator 

of the scalar field and the role of the pion loops in this propagation : 

2 



nuclear 

Xs 



-2 



2 

vac 



1 



a 



(tt) 
N 



+ CT 



N 



J 77 



*2 



m 



*2 



cr 



N 



n5s(o) 



1 



m 



*2 



(51) 



Here 1155(0) is the full scalar polarization propagators (in which we include the coupling 
constant) : 

M* 

^ssiO) = gf no(0) [1 - Vres TiMV' ■ , (52) 

Ep 

where Kps is the residual interaction : 



i2^FC 



+ 



M*N 9s 



*2 



E*p mf 



(53) 
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This expression differs from the mean-field one of our previous work [5] by the presence 
of the term in d'^e^^ / dp^ . In the expression [51] we have introduced quantities that we 
denote effective sigma commutators with the exphcit expressions : 



2\FC /!#* ^2 



ON)en = {a';%ff + {a'^%ff = cr^^' + '^^-^^ + "^2- (54) 

It is interesting to compare it with the effective sigma commutator governing the evolution 
of the chiral condensate. which is : 

[(^Nhot = {(Tn) + {(^n) = —^n + -7. — + /-"^^-- (55) 

P ^ P P 

As was shown in ref. [7j, this last quantity, corresponding to the full nuclear sigma 
commutator per nucleon, also governs the evolution of the pseudoscalar susceptibility 
according to : 

,,,= _2M%fl- (Mi^V (56) 



<> 7T 7T \ <> TV 7T 

Comparing the expressions [51] and [55] we see that they are not identical. They coincide 
to leading order in density but the higher order terms differ. In particular ($^)'^*-^/p 
is replaced in eq. [54] by d{^'^)^'" /dp. We recover in this work the result of [8J, that 
the quantity which governs the transformation of the quark scalar density fluctuations 
into nucleonic ones represents an effective nucleon sigma commutator which includes pion 
loops. In ref. [Bj it was the free nucleon one while here interactions are incorporated and 
we deal with effective values. Our present result thus generalizes the first order result of 

m 

Coming back to the expression [5l] of Xs"^'*^"^? we point out that the pion loops only 
enter in the effective sigma commutator which includes medium effects. The mass, m*, 
which enters the denominators is instead unaffected by the pion loops. The nuclear 
contribution ^^g^^/ear indeed linked to the propagation of the chiral invariant scalar field 
s which is not coupled to pions. The pure pion cloud contribution (eq. [H]) is totally 
decoupled. This point was already discussed, within a different approach in ref. [12] . 



4 Numerical results and discussion 



4.1 Fixing the parameters 

One important constraint to fix the parameters is the pion-nucleon sigma term which 
has a value, aN — 50MeV. In our approach it receives two contributions: one from the 
pion cloud, the other from the sigma meson. We point out that if our model represents 
a bosonized underlying NJL model, the latter contribution represents the sigma term of 
constituent quarks. It is interesting to make a connection with the lattice evaluations 
of the evolution of the nucleon mass with the pion mass. Indeed lattice simulations of 
the nucleon mass as a function of the squared pion mass (equivalently the quark mass) 
are available in the mass region beyond m,r — 400MeK. The derivative dMjq / dmi. = 
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(TAr/m^ provides the nucleoli sigma commutator. In turn the derivative of leads to 
the susceptibility. Both quantities are strongly influenced by the pion cloud which has 
a non-analytic behavior in the quark mass, preventing a polynomial expansion in this 
quantity. However, it is possible to extrapolate lattice data using chiral models of the 
nucleon as discussed by Thomas and collaborators and we use their recent version ref. 
[2T] . In their work the pionic self-energy contribution to the nucleon mass is separated 
out using different cut-off forms for the pion loops (gaussian, dipole, monopole) with an 
adjustable parameter A. They expanded the remaining part in terms of m\ as follows : 

MAr(m^) = ao + a2ml + a^ml + ^^^(m^. A) . (57) 

The best fit value of the parameter which fixes the susceptibility shows little sensitivity 
to the shape of the form factor, with a value 04 ~ —0.5 GeV~^ while 02 — 1.5 GeV~^ (see 
ref. [21])), from which we can infer the non-pionic pieces of the sigma commutator : 

^non-p^on ^ ^2 ^ ^ ^2 ^ 3 04 ~ 29 MeV . (58) 

This number indicates the existence of a large component beside the pion one, that 
it is natural to attribute to the scalar field. The identification of 0-^°""^*"" with a^'' of 
our model fixes the sigma mass to a value rua = 800 MeV, close to the one ~ 750MeV 
that we have used in our previous article [5]. As it is the ratio gs/fn^ which is thus 
determined this value of rricr is associated with the coupling constant of the linear sigma 
model gs = MN/fn = 10. Lowering gs reduces nia. With our cutoff , A = 0.98 Gel^, 
which yields cr^'' = 21.5MeV, the total value of the sigma term is = 50.5 Mel^, a 
quite satisfactory result. 

Although the procedule becomes more uncertain, on can also try to extract the non- 
pionic part of the nucleon scalar susceptibility : 

non— pion r,\^^/vac ^ I \ \HHlvac a lKCt\ 

<> TV ^\ / 

It has to be compared with the non-pionic piece of the nucleon scalar susceptibility. This 
can be done taking the zero density limit of eq. [51] : 

xK = -2^ (4i - ^) - = -2^T - ■ m 

This identification allows an estimate of the crucial parameter k^vs- For convenience, as 
in [5], we introduce the dimensionless parameter C = {f^/2Mj^)Kiys and we obtain : 

a, = -^j3-2C). (61) 

The value 04 ~ —0.5 GeV~^ leads to C ~ 1.25, which indicates a strong cancellation 
between the two components of 04, which makes the results sensitive to the exact values 
of the parameters. Therefore, due to the uncertainties of the procedure, we take this value 
of C only as indicative. 
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UJ 




Figure 1: Binding energy of nuclear matter with g^^ = 7.3 and C = 0.985 keeping only 
the Fock term on top of a and uj exchange. The full line corresponds to the full result, 
the dotted line represents the binding energy without the Fock term and the dot-dashed 
line corresponds to the contribution of the Fock term. 



Hence for the following discussions we fix the various parameters according to rrirj = 
800 Mel^, A = 980 Mel/, c/^^ = 0.7, g'j^^ = 0.3, g'^^ = 0.5, Cp = 2 and the parameter 
associated with the nucleonic response is allowed to vary in a broad region around C = 
1.25. The vector coupling constant g^^ is a totally free parameter. Finally, as in ref. [5j, 
we also impose that the nucleonic response vanishes at the chiral restoration [s = —fj^) : 



kns{s) 



N 



ds 



o2 



S 



Accordingly, the effective nucleon mass of eq. [3] is modified according to : 



M 



N 



1 



and the effective scalar coupling constant becomes dependent on s: 



1 



ds 



+ I^NSS 1 



2/. 



(62) 



(63) 



(64) 



4.2 Results with the influence of the Fock term 

The result of the calculation for the binding energy of nuclear matter in presence of the 
Fock term only (i.e., ignoring the correlation energy) is displayed in fig. [H The saturation 
point(po = 0.16 /m-3, E/A = -16.1 Mel/) is obtained with C = 0.985 and g^ = 7.3. 
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Figure 2: Density evolution of the nucleon (dot-dashed curve) and sigma masses (full 
curve) with g^^ = 7.3 and C = 0.985; the dashed curve corresponds to the sigma mass 
when the effect of the nucleon susceptibility is removed. 
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Figure 3: Density evolution of the QCD susceptibilities (normalized to the vacuum value 
of the pseudoscalar one) with = 7.3 and C = 0.985 keeping only the Fock term on 
top of (7 and u exchange. Dashed curve: pseudoscalar susceptibility. Full curve: scalar 
susceptibility. Dotted curve: nuclear contribution to the scalar susceptibility. Dot-dashed 
curve: pion loop contribution to the scalar susceptibility. 



16 



60 



50 >: 



40 - 



CO 

5 

S 30 



(f> 20 



10 



I ' ' ' ' ' 1 

0.5 1 1.5 2 2.5 3 

P/Po 

Figure 4: Density evolution of the full in-medium sigma term, {0^)101 (dotted line) and 
of the effective sigma term, {crN)eff (full line), with = 7.3 and C = 0.985 keeping only 
the Fock term on top of a and u exchange. 
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Figure 5: Density evolution of the residual interaction in MeV'"^ (continuous line). 
Dashed line: omega exchange contribution. Dot-dashed line: Fock term contribution. 
Dotted line: sigma exchange contribution. 
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The corresponding incompressibility is K = 254 Me V^. Comparing with our previous 
results [5] we see that the inclusion of the Fock term somewhat improves the description 
since it lowers K. We observe in fig. [T] that the value of the Fock term at saturation 
density is nearly the same as the binding energy. It turns out that, at this density, this 
binding comes from the transverse Fock term [g' + rho exchange). In the longitudinal 
channel there is an almost exact compensation between the pion and the g' exchanges. 
The genuine Fock term associated with the Yukawa piece of the pion exchange represents 
an attractive contribution to E/ A of only —4.2MeV ~ 20% smaller than what is obtained 
in relativistic theories ^24j where the irNN form factor is omitted. 

For completeness we also display in this subsection (fig. [2]) the nucleon and sigma 
masses with these values of the parameters, although they are not sensitive to the inclusion 
of the Fock term. At normal nuclear matter density the effective nucleon mass is close 
to 800 MeV. As already emphasized in our previous work [5] the sigma mass remains 
instead quite stable, due to the introduction of the nucleonic scalar response. 

The evolution of the susceptibilities is depicted in fig. [31 The convergence effect 
between the pseudoscalar susceptibility (dashed line) and the scalar one (full curve) is 
more pronounced than in the mean field approach of our previous work [5] . The reason is 
clear and linked to the introduction of pion loops. Firstly the evolution of the pseudoscalar 
susceptibility, Xps, as written in eq. [56l is governed by the in-medium sigma term, {aN)tot, 
which in the present approach receives a contribution from the pion loops which increases 
its value (see eq. [55]). Secondly, for the scalar susceptibility, the effective sigma term, 
{cN)eff, transforming the scalar quark density fluctuations into nuclear excitations also 
has a pion cloud contribution (eq. and this effect increases the nuclear part, ■^'^^'^'■^"■^ 
(dotted curve), of the scalar susceptibility. Finally the scalar susceptibility also contains a 
pure pion loop contribution, y^P^°^^°°P (dot-dashed curve), which becomes more important 
at large density. For illustration we also compare in flg. [Hthe in-medium modifled nucleon 
sigma term {aN)tot governing xps (dotted curve) and the effective sigma term {crN)eff (full 
curve) affecting the nuclear piece T^g'^'^'^"'". 

The density behavior of the residual interaction (eq. [531) is displayed in flg. [51 with 
separate curves for the various components (omega and sigma exchanges, Fock contribu- 
tion). While the u exchange has a very smooth density dependence, one notices the rapid 
variation of the sigma component associated with the dropping of the effective scalar 
coupling constant g^ with increasing density. Around Pq Vres turns from attraction into 
repulsion. This feature comes in part from the Fock term evolution and mostly from the 
behavior of the sigma exchange contribution 

4.3 Influence of the correlation term 

The expression of the longitudinal contribution to the correlation energy, which is always 
attractive, has been given in eq. [201 It is is depicted schematically in flg. [6land involves the 
full longitudinal spin-isospin polarization propagator 11^, which is calculated in practice 
within the RPA ring approximation. There is a similar expression for the transverse piece. 

With the previous choice of these parameters it turns out that the numerical value 
of the correlation energy is close to the Fock term energy. This extra attraction has to 
be compensated by an additionnal repulsion in order to account for the saturation point. 
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Figure 6: Schematic representation of the longitudinal spin-isospin contribution to the 
correlation energy. 
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Figure 7: Binding energy of nuclear matter with g^^ = 8, rricr = 850 MeV and C = 
0.985 with the Fock and correlation energies on top of a and uj exchange. The full line 
corresponds to the full result, the dotted line represents the binding energy without the 
Fock and correlation energies and the dot-dashed line corresponds to the contribution of 
the Fock terms. The decreasing dotted line (always negative) represents the correlation 
energy. 
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Figure 8: Density evolution of the longitudinal (full line) and transverse (dashed line) 
contributions to the correlation energy for g'^^ = 0.3. 
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Figure 10: Density evolution of the QCD susceptibilities (normalized to the vacuum 
value of the pseudoscalar one) with = 8, rricr = 850 MeV and C = 0.985 with the 
Fock and correlation terms on top of cr and u exchange. Dashed curve: pseudoscalar 
susceptibility. Full curve: scalar susceptibility. Dotted curve: nuclear contribution to the 
scalar susceptibility. Dot-dashed curve: pion loop contribution to the scalar susceptibihty. 
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Figure 11: Density evolution of the full in medium sigma term, {aN)tot (dotted line) and 
of the effective sigma term, {(JN)eff (full line), with g^^ = 8, mo- = 850 MeV and C = 0.985 
with the Fock and correlation terms on top of a and uj exchange. 
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Increasing only the vector coupling constant tends to give too low a saturation density. 
However a slight adjustment of the sigma mass yields a correct saturation point. Taking 
g^j = 8 and = 850 MeV (corresponding to = 26 MeV) and keeping C = 0.985 

we obtain the saturation point : po = 0.16 /m~^, E/A = —15.7 MeV . The corresponding 
correlation energy is E'-''°^^/A = —17. 4 MeV which divides into a longitudinal piece, 
E^orr^^ = -7.9 MeV and a transverse one, £;^°"7A = -9.5 MeV. The resuhs of the 
calculation are shown in fig. [71 We stress the importance of the transverse channel in 
the correlation energy as compared to the longitudinal one. The reason is the strong 
screening of the pion exchange by the short-range interaction in the longitudinal channel, 
hence a sensitivity to the g' parameters. For instance changing g'^^ from 0.3 to 0.5 reduces 
further the longitudinal value to E^"^^/A = —3.7 MeV while the transverse contribution 
shows less sensitivity, with E^"^^/A = —7.8 MeV , such that the total value becomes 
ECorr^j^ = —11.5 MeV. With the value 5?^^ = 0.3 the longitudinal contribution to the 
correlation energy is dominated by 27r exchange with one Delta intermediate excitation. 
Iterated pion exchange is smaller, because of the screening effect. Therefore, as the main 
mechanism with one Delta is a part of the A^A^ potential, we expect the longitudinal 
correlation energy to be linear in density, at variance with the transverse channel where 
the Delta is not dominant. These behaviors are illustrated in fig. [HI The longitudinal 
correlation energy has indeed a linear behavior at low densities. It deviates from linearity 
at high density where it becomes more attractive. This reflects the appearance of many 
body forces due to the increase of the nuclear pionic field, the critical opalescence effect 
[26] induced by the small value of (^^a- Taking for illustration a larger g' value, g'j^^ = 0.5, 
the linearity is fulfilled in a larger density range, as shown in fig. [9] 

It is interesting to make a comparison with the chiral perturbation calculation of ref. [9] 
where the effect of short-range correlation is absent. In the latter work the contribution of 
the iterated pion exchange is found to be —68 MeV, the corresponding Hartree diagramm 
coinciding with the lowest order term of the longitudinal correlation energy depicted on 
fig. [6l There is nearly one order of magnitude difference between the two results, owing 
to the fact that in our phenomenological approach pion exchange is strongly screened by 
short-range correlations. 

Although in detail our numerical results on the correlation energy present a sensitivity 
to various factors, the values of the g' parameters, the rho coupling (Cp) or the form 
factors, especially in the transverse channel, we believe that our conclusion about the 
moderate importance of the correlation part of the energy is robust. The important point 
is that the rest of the interaction, as provided by lattice data, is sufficient to bind the 
nuclear system. 

In the spirit our approach has a similarity with the density functional approach of 
Finelli et al [16]. As them we have background scalar and vector fields. The first one 
has a connection with chiral symmetry restoration of QCD. However there is an essential 
difference concerning the scalar field. For Finelli et al. the full nucleon sigma commutator 
enters the scalar self-energy of the nucleon according to : 

(0) CTtvMat 

In our case instead, only part of the nucleon sigma commutator enters this self-energy. It 



22 



is the part which originates from the chiral invariant field s. To leading order in density 
we have : 

= m' = -^-p^ PS. (66) 

The other component from the nucleon pion cloud, a^^'^ ^21 MeV should not enter the 
mass evolution as it is cancelled by other terms, as imposed by chiral constraints. This 
point was emphasized in several works, see e.g. reference [271 fT2]. Accordingly the link 
between our nuclear scalar field and QCD does not occur through the total nucleon sigma 
commutator but only through part of it, its non-pionic piece. Information on this last 
quantity can be obtained from the lattice result on the nucleon evolution with pion mass, 
once the pionic contribution is separated out, as done in ref. [21] and explained in section 
4. Our result is therefore not totally model independent but this limitation is imposed 
by chiral constraints [28] . 

For the evaluation of the scalar susceptibility including correlation terms, we make the 
assumption that the pionic correlation energy density and the correlated part of the pion 
scalar density go both like p^, a reasonable approximation in view of the results shown 
on fig. [HI The net conclusion is that the convergence effect between the scalar and pseu- 
doscalar susceptibilities is much more pronounced as shown in fig. [10] (compare with fig. 
[3]). One reason is that the pionic piece of the scalar susceptibihty, ■y^P^°^^°°P ^ jg significantly 
increased by the RPA correlations, a conclusion already reached in our previous work 
[T2] . Another reason comes from the fact that the full in-medium sigma term {d'i\f)tot, 
governing the pseudoscalar susceptibility and the effective sigma term, {(JN)eff affecting 
the nuclear piece of the scalar susceptibility remain relatively stable with the density, as 
depicted on fig. [Til while they significantly decreased in figJU 



5 Conclusion 

We have studied in this work a relativistic nuclear model based on a chiral version of 
the a and u exchange model. We have worked beyond the mean field approximation 
introducing the pion loops. Our goal is to reach a consistent description of matter which 
can apply at densities larger than the saturation one by introducing the pion loops on 
top of a mean field approach. The aim is to reach the best possible description in the 
framework of the model. The pion is an important actor of the nuclear dynamics and a 
credible theory should incorporate its effect. 

In establishing the model and fixing its parameters we have kept contact with other 
domains. In particular we introduce beside the pion the short-range components, which 
are indissociable from the pion in the spin-isospin residual interaction. They are em- 
bedded in the Landau-Migdal parameters and for their values we have used the most 
recent informations. They favor a clear deviation from universality, with a small value of 
g'^^ ~ 0.3. We have also insured the compatibility of our parameters with informations 
from QCD. One important constraint is the total nucleon sigma commutator. But more 
information is also available from the lattice results on the evolution of the nucleon mass 
with the quark mass. From the analysis of these data by Thomas et al. [21] we extracted 
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the non-pionic piece of the sigma commutator. In our model this is the scalar field contri- 
bution and it fixes the sigma mass to be used in our inputs. Another piece of information 
is obtained from the higher term of the expansion, even if wc take this information as 
indicative. Combining the informations from spin-isospin physics and from QCD lattice 
results we can reach a tenable description of the saturation properties. /^From spin-isospin 
physics, we have taken the values of the Landau- Migdal g' parameters to be used in as- 
sociation with pion and rho exchanges. They limit the value of the correlation energy 
from the pion loops through a suppression of pion exchange by the short-range compo- 
nent. On the other hand the value of the non-pionic piece of the sigma commutator from 
lattice QCD, which fixes the sigma mass, provides the attraction needed for saturation 
to occur. Moreover the higher term 04 of the nucleon mass expansion strongly limits the 
many-body effects that can occur in the propagation of the nuclear scalar field. Other- 
wise the tadpole terms chiral theories alone would induce too much attraction and destroy 
saturation. It is quite remarkable than the parameters which reproduces the saturation 
properties are compatible to the lattice ones. Nowhere have we faced a contradiction. It 
was not a priori obvious that the amount of attraction needed from the scalar field fits 
the value of the non-pionic sigma commutator from lattice data. Nor that the amount of 
cancellation necessary in the medium effects of the sigma propagation from the nucleon 
scalar response is reproduced in the expansion of the lattice results, with values of the 
parameter C remarkably close (C = 1.25 for the lattice results versus C ~ 1 in our fit). 
This consistency gives credibihty to the description.. 

In more details, we have investigated the role of the Fock term and of the correlations. 
For the second one we have found a moderate value, ~ — 17 MeV at po due to a suppression 
of the pion contribution by the short-range component in such a way that the transverse 
channel linked to rho exchange becomes dominant. The longidudinal component of the 
correlation energy is linear in density but it evolves faster at high densities, signalling 
the possible appearance of many-body forces of attractive nature. They arise from the 
enhancement of the nuclear pion field, the critical opalescence efi^ect which softens the 
equation of state at large densities. The compatibility of the parameters of our approach 
with the QCD lattice data is a support for our results on the density evolution of the 
QCD quantities, the quark condensate and the QCD suceptibihties.We have insured that 
these results are in addition fully compatible with the saturation properties of nuclear 
matter. For the quark condensate we find small deviations from a linear behavior with 
density, in spite of the inclusion of various types of many-body effects. This also applies 
for the pseudoscalar susceptibility which follows the condensate. The scalar susceptibility 
is instead very sensitive to the interactions and shows a large enhancement as compared 
to the free value. It even surpasses the pseudoscalar one beyond p ~ 2.2 po while in the 
vacuum it is much smaller. 

Our description could be improved in various directions such as the incorporation of all 
the Fock terms (sigma and omega exchange, time component of the rho meson exchange) 
and a covariant description of the short-range interaction in a fully covariant framework. 
It would also be interesting to extend the approach to the case of asymmetric nuclear 
matter and neutron matter. 
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